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We derive threshold equations for self-organization of laser driven atoms in an optical cavity. Our 
analysis includes probing with either a traveling wave or a retro reflected lattice. These two scenarios 
lead to qualitatively different behavior in terms of the response of the system as a function of cavity 
detuning with respect to the probe. In addition our analysis includes the effects of an intra-cavity 
trapping potential which is also shown to impact on the threshold condition. We specifically consider 
the case of an intra-cavity lattice as used in [1[ but our treatment can easily be modified to other 
geometries. 



I. INTRODUCTION 

Atoms coupled to the standing wave mode of a cav- 
ity will self-organize when driven by a sufficiently intense 
transverse laser field [2j . As the atoms scatter more light 
into the cavity, a potential is produced which localizes 
the atoms in a more favorable configuration for scatter- 
ing. This further enhances the localization of the atoms. 
Thus above a threshold pump intensity, an initially uni- 
form distribution of atoms will undergo a phase transi- 
tion, spontaneously reorganizing into a lattice configura- 
tion. Equations determining the threshold intensity are 
given in Q but they are specific to a set up in which 
the atoms are driven by a standing wave field. They 
are also restricted to a particular probe detuning, and 
an initial distribution of atoms in which its transverse 
extent can be neglected. Here we extend the analysis 
to include the effects of an external confinement of the 
atoms and a traveling wave probing configuration. We 
restrict our attention to a lattice confinement as used 
in recent experiments [l| but the analysis can easily be 
modified to other trapping configurations. For the case 
of a traveling wave probe it is found that the threshold 
behavior is qualitatively different to that of a retrore- 
flected lattice geometry and is analogous to a collective 
atom recoil laser (CARL) [4[. We first give an overview 
of the self-organization process introducing the needed 
terms and equations for our analysis. We then give a de- 
tailed derivation of the threshold conditions for both the 
traveling wave and retro-reflected probing conditions. 



II. SELF-ORGANIZATION OF THERMAL 
ATOMS IN EXTERNAL POTENTIALS 

To derive the equations for the onset of self- 
organization of atoms coupled to an optical cavity, we 
extend the previous theoretical treatment of Q to include 
an external trapping potential and the spatial extent of 
the atomic distribution transverse to the cavity. In ther- 
mal equilibrium the spatial density distribution of atoms 



in a confining potential V(x) is given by 

*w-i«p(-$), <u 

where Z = j exp (— V(x))/(fcBT)) d 3 a;. For our case the 
confining potential consists of two components 

V(x) = ^r(x) + V s (x), (2) 

where Vt(x) is the potential associated with an exter- 
nal potential and Vs(x) is the potential associated with 
the probe beam itself and light scattered from the probe 
beam into the cavity. 

Self-organization in this system arises from the back- 
action of the scattered light on the atomic distribution. 
As the atoms scatter light into the cavity, the resulting 
potential Vs(x) localizes the atoms to a more favorable 
configuration for scattering. This results in increased 
scattering which then further enhances the localization 
of the atoms, and, under sufficient driving strength, an 
initially uniform distribution of atoms will undergo a 
phase transition spontaneously reorganizing into one of 
two possible lattice configurations. 

In our analysis of the problem two spatial averages play 
an important role, namely 

a = J f(r,z)p a (x)d 3 x (3a) 

and 

/3 = J f(r,z)h(x)p a (x)d 3 x, (3b) 

where, f(r, z) — cos(fcz) exp(— r 2 /w 2 ) is the mode func- 
tion of the cavity and h(x) the mode function of the probe 
beam, which in our experiments has either a standing or 
traveling wave form. The factor a accounts for the re- 
duced coupling of the ensemble to the cavity due the spa- 
tial positions of the individual atoms relative to the cav- 
ity mode. The factor (3 quantifies the scattering into cav- 
ity due to the spatial organization of the atoms. It takes 
the role of an order parameter having a value near zero 
prior to the phase transition and rapidly approaches ±1 
as the probe coupling, Q, rises above a threshold value. 



2 



Since the potential Vs in general is a function of a and 
/3, Eq [1] and [3] must be solved self-consistently. This can 
be done by first expressing Vs in terms of arbitrary a 
and {3 and using EqQ]to determines the atomic distribu- 
tion. The distribution function is then substituted into 
Eq [3] to give a coupled set of self-consistency equations 
for a and f3. To determine the threshold condition for 
self-organization to occur we note that, near the critical 
point, j3 ~ 0. Hence we can expand the self-consistency 
equations to lowest non-zero order in /3. This yields the 
necessary condition for there to be a non-zero value of /3. 

We first determine Vg (x) by looking at the steady state 
solutions describing the system for arbitrary a and /3. At 
large atomic detunings, we can adiabatically eliminate 
the atomic excited states. The resulting Hamiltonian for 
a collection of two- level atoms with a fixed density, p a (x), 
inside a cavity driven from the side by a classical field is 

H = (-A c + NU a)a)a + N(Sl R /3)*a + N(n R p)a) (4) 

with the cavity decay described by the Liouvillian 

Cp = (k + NToa)(2apa ! — a/ap — pa) a). (5) 

In these expressions A c = uj p — uj c is the cavity detuning 
with respect to the probe, Uq is the single atom disper- 
sive shift, il R is the pump rate per atom of the cavity due 
to scattering of the probe field, and To is a correction to 
the cavity decay associated with spontaneous emission. 
In terms of the atom-cavity coupling, g, atom-probe cou- 
pling, f2, the detuning, A, of the probe from the atomic 
resonance, and the atomic linewidth, 7, the parameters 
Uq, To, and Q, R may be written 



In our system ATo <C k and so we neglect this term in 
all that follows. 

In steady state the cavity mode is described by a co- 
herent state with amplitude 
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(6) 



where A c = A c — NU^a is the detuning from the disper- 
sively shifted cavity. 

The total scattered potential, Us(x), is then given by 



Using the resulting potential V(x) = Vr(x) + Vs(x) to 
determine the density distribution via Eq. [1] will lead to 
self-consistency equations for a and (3 via their definitions 
in Eq. |3] 



III. DERIVATION OF THE THRESHOLD 
EQUATIONS 



by 



The self-consistency equations for a and /3 are given 



a = 1 f e -(^T(x) + Vs(x))/(fe B T) /2(? , z) d 3 x (Qa) 



and 



13 = ± I e -(V T (x)+y S (x))/(fe B T) /(r) z)h{x) d 3 x (gb) 



with 



Z= / e -(VT(x)+Vs(x))/(fc B T) d 3 x _ 



(10) 



We henceforth consider the specific case in which Vt is 
given by an intra-cavity FORT that has a wavelength 
twice that of the probe as used in recent experiments [l| . 
This traps the atoms at every alternate antinode of the 
cavity . In this case, the integration along the z direction 
can be carried out by integrating over single trapping 
sites FORT lattice potential and weighting each site in 
accordance with the probability that the site is occupied. 
Since the integrands have the same periodicity as the 
lattice potential, the integrations over each trapping site 
are equal and we need only consider a single site. This 
is equivalent to considering all the atoms to be located 
in one site. Furthermore, we assume the atoms to be 
sufficiently well localized by the FORT potential so that 
the trapping site can be well represented by its harmonic 
approximation. We thus replace the potential Vt by its 
harmonic approximation 

r r ^ 22 22 

Vt — -muj r r + -mu),z 
2 2 

and extend the z integration out to infinity. 



A. Lattice geometry 



V s (x) = -\nh(x) + \gf(x)\ 2 



where 



\h(x)\ 2 + 2e/(r, z)Re (/Te^ 9 h(x)) 



+ e 2 \f3\ 2 f(r,z) 

NUn 



and 



A? 



A? 



(7) 



(8) 



In the case of probing the system with a lattice, the 
mode function is 



h(x) = cos(kx 



(11) 



where k = 2tt/X and A the wavelength of the probe light. 
The potential V s is then given by 

nn 2 



cos (kx 



+2e/3/(r, z) cos(kx - <p) cos(6») 



(12) 



3 



where e and are determined by Eq [8j Note that we 
have dropped a term proportional to /3 2 as we will only 
be interested in terms up to first order. When we expand 
Eq[5]to lowest non-zero order in (3 we encounter integrals 
of the form 



2 = 



-* 2 /(2- 2 ) 



W(cos(kx — (j)))dx 



for some function W . Here a is approximately the spatial 
extent of the atomic distribution along x. These integrals 
can be greatly simplified when ka 3> 1 which means the 
spatial extent of the atomic distribution is larger than the 
wavelength. In this case the integrand can be averaged 
over the wavelength. Denoting the averaging by () we 
then have 



(W(cos(fcz; - 0))) 



(W{cos(kx))) crV27r 



(13) 



Thus, within this approximation, the phase </> of the lat- 
tice potential is irrelevant and we subsequently set it to 
zero. The expansions of Eq[9]can now be readily evalu- 
ated. 

The zeroth order expansion of Eq [9a| gives 



a 



where 



J e -V T I {k B T) ei i cos 2 (fcr) z j £3 X 

J g—Vr /(kBT) g/i cos 2 (kx) ^3^. 

/e-^/(^ T >; 2 (r,z)d 3 x 
J e -v T /(k B T) d 3 x 



fi = ~nn 2 /(Ak B T) 



(14) 



is the depth of the potential due to the probe relative 
to the thermal energy k B T of the atoms. Within the 
harmonic approximation for Vt(x) we then find 
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-4/r, 



2/rj 



(15) 



where 77 = V T a/(k B T) and V T0 is the depth of the FORT 
potential. Note we have utilized the fact that the waist 
of the FORT beam is V2 larger than that of the cavity 
mode associated with the probe wavelength. 

For the expansion of Eq I9b[ the zeroth order term is 
proportional to an integral of the form 



-x 2 / (2a 2 ) u cos 2 (kx) 



cos(fcx) dx. 



Within the approximation ka 3> 1 this term is zero as 
expected for nearly uniform distribution per wavelength. 
If the approximation ka ^ 1 is not satisfied then the 
zeroth order term would be non-zero and /3 would be non- 
zero for any value of the probe coupling. This is because 
the atoms would already have sufficient localization to 



provide some level of scattering into the cavity for any 
probe intensity and thus no threshold would exist. 
The expansion of Eq [9b] to first order then gives 

P = 2efj,/3cos(6) 

J e -V T /(k B T) e „co S 2 (kx) cos 2(kx)f( r , z) & z x 



2efif3 cos(9)a 



J g — Vt I (k B T) g/i cos 2 (kx) ^3 j. 
^cos 2 (kx) cos 2( fca ,)^ 



e^cos(6')a 1 + 



^g/i cos 2 (kx) ^ 

Wg)Y 







(16) 



A nontrivial solution for f3 requires the expression within 
the square brackets to be unity thus giving the required 
threshold equation 
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NU a A, 



(17) 



It is worthwhile comparing our result to the original 
treatment in [3j . Our threshold equation is derived under 
the assumption that the atoms are trapped at every sec- 
ond antinode of the cavity and have a transverse spatial 
extent that is large relative to the wavelength. The result 
given in Q on the other hand neglects the transverse di- 
mension and assumes a uniform distribution along the 
z-axis. However these differences do not significantly im- 
pact on the threshold equation. Confining the atoms to 
every second antinode merely alters the parameter a and 
the net effect is a slight rescaling of the threshold value. 
Including the spatial distribution of the atoms transverse 
to the cavity gives rise to the ratio Ii(fi/2)/ Iq{^l/2) in the 
threshold equation. This ratio should be set to « 1 if the 
atoms are confined to length scales much smaller than 
the wavelength along this dimension. Thus recovering Eq 
(19) of ref Q amounts to setting a = 1/2 to account for 
the initially uniform distribution of atoms along the cav- 
ity axis, and setting 7i(/x/2) / 'Iq(^l/2) = 1 to account for 
neglecting the transverse dimension. Using their chosen 
detuning of A c = — k then gives exactly the expression 
reported in ref 0]. Our threshold equation is therefore a 
generalization of the result given in [3j to account for an 
arbitrary detuning of the cavity and for the transverse 
extent of the atomic distribution. 



B. Traveling wave geometry 

When probing with a travelling wave beam, the mode 
function is 



h(x) = 



— ikx 



(18) 



and the potential V s is given by 



V„ = 



hn 2 



2e/(r,z)Re(/3V ( ' ca; - e) ) 
[2e/(r,z)|/3|cos(fcc-0)] 



(19) 
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where 6 = 9 + arg(/3). Note that we have dropped the 
constant term hft 2 /A and as before we have neglected 
the term proportional to /3 2 as we are only interested in 
terms up to first order. 

Taking the zeroth order expansion of Eq [9a] we get ex- 
actly the same expression for a as we do for the lattice 
case. This is not surprising as a represents the averaging 
of the dispersive shift due to the spatial distribution of 
the atoms. As we assume the spatial extent in the trans- 
verse direction is much larger than the wavelength this 
averaging will be unaffected even if a lattice potential is 
invoked along this dimension. 

For the same reason as for the lattice case, the zeroth 
order term in the expansion of Eq|9b] vanishes and to first 
order we have 



P = 2e/*|/9| 



/ e-^/( fcsT )/ 2 (r, z)e tkx cos(fcx - 



i d 3 x 



= 2en\l3\a(e l 



J e -v T /(k B T) d 3 x 
: cos(kx -0)) 



e[i\p\ae 
(e^ae 



(20) 



Equating the term in parentheses to unity then yields the 
equation 



A 4 = 



A? 



-NUqol 



(21) 



However, since /i is real, this equation has no solution due 
to the e %e term. The difficulty stems from the fact that, 
for an organized array of atoms in the cavity, the scat- 
tered potential has minima that are displaced from the 



positions of the atoms. The result would be that, as the 
transverse lattice potential forms, the atoms would expe- 
rience a force along the direction of the probe, establish- 
ing a moving lattice similar to the situation encountered 
in the collective atom recoil laser (CARL). We therefore 
argue that the appropriate threshold equation for this 
geometry is Eq [5T] with the e 10 term omitted. We also 
note if an appropriate force were applied to the atoms, 
they would be held at a position offset from the poten- 
tial minima giving a self consistent solution for a and /3. 
Such a force could be provided by the Gaussian profile of 
the cavity mode and would have the effect of stabilizing 
the lattice by canceling the phase term in Eq [2T1 



IV. CONCLUSION 

We have derived threshold conditions governing the 
onset of self-organization for N atoms in a high finesse 
cavity. Our equations generalize those given in Q by al- 
lowing for arbitrary detuning of the probe relative to the 
cavity, and account for the full spatial extent of the initial 
atomic distribution. We have also included the case of a 
traveling wave probe that has not yet been considered. 
In this case the response of the system is qualitatively 
different from that of a lattice probe. Specifically, self- 
organization occurs for a probe detuned above the cavity 
resonance contrary to the lattice case. Our derivations 
have been restricted to the specific case which includes 
an intra-cavity FORT but the treatment can be easily 
modified to other external trapping potentials. 
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